Teaching the Computer how to Discover(!) and then Prove(!!) (all by Itself(!!!)) 
Analogs of Collatz's Notorious 3x+l Conjecture 



Doron ZEILBERGER^ 
Mathematics: an Experimental Science 

In spite of their exponential growth, computer-assisted and especially computer-generated mathe- 
' matical research are still in their infancy. One approach is that of formal proofs using the axiomatic 

' method. This method is based on the myth, that goes back to Euclid, that mathematics is a de- 

ductive science, where one starts with a bunch of axioms, (initially supposed to be "self-evident", 
but later conceded as true-by-fiat) and then uses rules of deduction, and step-by-step, arrives at 
" . (seemingly) non-trivial results. 

^ ; 

' Of course, mathematics could be presented that way, and unfortunately, often is. But that is not 

O""""' how it is discovered. Pretending that the axiomatic method is how mathematics should be done, 

. and trying to indoctrinate poor computers to do it that way, is a highly inefficient use of computers' 

• . time. 

^ ' Deep inside, mathematics is, or at least should be, an inductive science. How do human mathe- 
, ^ maticians come up with such amazing conjectures? By experimenting! How do they come up with 
such amazing proofs? By experimenting!, and not by combining axioms. 

> • 

' For many results in mathematics, one only needs one non-trivial "axiom" , that of Peano's "axiom" 

' of induction: 

O ■ 

; P(0) & (P(n) ^ P{n + 1)) ^ P{n) for all n > . 

cn ; 

, All we have to do then is to have the computer prove (all by itself) -P(O) (an a priori routine fact), 

i and the not-a-priori routine P{n) =^ P(n + 1), but, with the help of symbolic computation, treating 

^1 n as a symbolic integer (which is one object) rather than as a variable (ranging over positive 

\^ ' integers), we might hope to have the computer do it for us. 

. . , Often one tries to prove a statement by induction, but fails. In that case one has to try again. If 

P{n) is inadequate to prove P{n + 1), perhaps we need another statement Q{n), such that Q{n) 
and P{n) imply P(n + 1). Alas, now we also have to prove that P{n) and Q{n) imply Q{n + 1). 
This may be easier-said-than-done, and we may be forced to introduce yet-another statement R{n). 
If we want to avoid a Ponzi scheme, we need this process to finally halt, but if we can train our 
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computer how to do it, it will be able to keep trying, before halting or giving up, for a much longer 
time. 



Human mathematicians also do symbol crunching, most of which is purely routine, and that is 
much better delegated to computers. 

Research is one percent Inspiration and ninety nine percent Perspiration 

At this time of writing, humans are still needed to find ideas and strategies for generating con- 
jectures, and for proving these conjectures. But once the human has some ideas, it is much more 
efficient to teach the computer these ideas, and let the computer search for conjectures, and most 
impressively, proofs, all by itself! 

In this case-study in computer-generated mathematics, I will describe how I read very carefully 
a beautiful mathematical paper [AGKL], written by four brilliant human beings: Ed Grove and 
Gerry Ladas, and their (at the time) respective students. Candy Kent and Amal Amleh. I then 
extracted the ideas, looked at the structure of the proofs, but ignored the details. I then taught 
(programmed) my beloved computer to execute these ideas, and since it is much more patient than 
a human being, it was able to prove many more results. A very preliminary output can be viewed 
at the webpage of this article 

littp : //www . math . rutgers . edu/"zeilberg/mamarim/mamar imhtml/ collatz . html , 

where there are 144 computer- generated results, that include those proved in [AGKL]. The computer 
did all the phases of mathematical research all by itself: conjecturing theorems, conjecturing proofs, 
and finally, verifying that the conjectured proofs are correct. 

Background: The 3x-|-l problem 

Mathematics abounds with easy-to-state yet (apparently) impossible-to-prove statements, for ex- 
ample the Goldbach and twin-prime conjectures, but none of them rivals the simplicity of statement, 
and most probably, difficulty of proof, of the Collatz infamous 3x+l conjecture, so beautifully ex- 
posited in Jeff Lagarias' [L] masterpiece. I am sure that all the readers know it, but for the sake of 
completeness, let me state it anyway. 



subject to the initial condition x(0) = xq. Then for every positive integer xq, the sequence is 
eventually periodic with the trivial cycle (2, 1). 

Paul Erdos claimed that mathematics is not yet ready for solving such problems. I strongly believe 



The Collatz 3x-|-l Problem 



Let x{n), (n > 0), be a sequence of positive integers defined by the first-order recurrence 




x{n) is even ; 
x(n) is odd . 
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that it very soon will be, thanks to the emergence of computer-assisted and computer-generated 

methods, that eventually would be able to rigorously prove such statements, for all xq < oo, as 
opposed to just numerically verifying it for .To < M for some large specific integer M. Of course, we 
can't do it naively, since there are infinitely many cases to check, and even our largest and fastest 
computers are finite (and so are we, and our universe, for that matter). 

If you try something and you fail, try and try again, then generalize and/or analogize], no use being 
a damn fool and trying to prove the original conjecture. In 1995, Clark and Lewis[CL] played with 
an analogous, second- order recurrence 



and started experimenting with it, numerically, by trying out some random initial values, and 
observing what is going on. For example, the following trajectory arises with x{—l) = 11, x{0) = 16: 



11, 16, 5, -11, -3, -7, -5, -6, -1, 5, 2, -3, -5, -4, 1, 5, 3, 4, 1, -3, -1, -2, -1, 1, 0, -1, -1, -1, -1, . . . 
while the following one, with x{—l) = 13, a;(0) = 6: 

13, 6, -7, -13, -10, 3, 13, 8, -5, -13, -9, -11, -10, 1, 11, 6, -5, -11, -8, 3, 11, 7, 9, 8, -1, -9, -5, -7, - 

1, 7, 4, -3, -7, -5, -6, -1, 5, 2, -3, -5, -4, 1, 5, 3,4, 1, -3, -1, -2, -1, 1, 0, -1, -1, -1, -1, -1, . . . 
or this one, with a;(— 1) = ll,a;(0) = 5 



They noticed that none of the trajectories they encountered seem to go to infinity, but instead 
eventually end-up either being identically 1, or identically —1, or end-up with the six-period 
(—2,1,3,2,-1,-3) for ever after. They then (probably) wrote a short computer program that 
checked their guess for all initial values —100 < x{—l),x{0) < 100, and they had ample empirical 
evidence for the following 

Theorem([CL]) Let x{—l) and a;(0) be any two integers whose largest odd common divisor is 1, 
then x{n) is either eventually the constant 1, the constant —1, or the six cycle (—2, 1, 3, 2, —1, —3). 

They then worked fairly hard to give a purely human proof. 

A few years later, in the remarkable paper [AGKL] already mentioned above, Amleh, Grove, Kent 
and Ladas considered all 16 recurrences, for a,P,j,S G {—1, 1}. 




x{n — L}-\-x{n) -r / 1 \ I / \ 

-^^ ^ — II x[n — 1) + x{n) IS even ; 

—x{n — l)+x{n), if x{n — 1) + x{n) is odd. 



11, 5, 8, 3, -5, -1, -3, -2, 1, 3, 2, -1, -3, -2, 1, 3, 2, -1, -3 . . . 




^-^ ^ if x[n — 1) + x[n) IS even; 

5x(n — 1) + 7x(n), if x{n — 1) + x{n) is odd. 



{AGKL) 



It is immediate to see that if the initial values x{—l),x{0) have an odd common factor, all 
the remaining terms would also, so without loss of generality, one can consider the case where 



x{—l),x{0) do not have a common odd divisor. Also if they have a common even divisor, as long 
as two consecutive terms are even, the transformation is "shrinking" , so without loss of generality 
gcd{xi-l),x{0)) = l. 

By a lot of human ingenuity and a considerable amount of human toil, Amleh et. al. proved that 
for two cases there are trajectories that go to infinity ((a, /?, 7, (5) = (1, 1, 1, 1), (1, —1, 1, —1)), and 
for ten other cases (out of the possible sixteen cases) they proved that all trajectories end-up with 
periodic orbits, and listed all of them. They left four open cases. One case is 

(a,/?,7,<5) = (-1,1,1,1) , 

for which they conjectured that all trajectories end up in one of the cycles 

(0,1,1) , (0,-1,-1) , (2,5,7,1,-3,-2,-5,-7,-1,3) , 

and the case 

(a, A 7, 5) = (-1,1,-1,-1) , 
for which they conjectured that all trajectories end up in one of the cycles 

(-1,-2,3) , (1,2,-3) , (-1,0,1,-1) , (1,0,-1,1) . 

Later they found yet another cycle, of length 33: 

(79, -31, -55, -12, 67, -55, -61, -3, 29, 16, -45, 29, 37, 4, -41, 37, 39, 1, -19, -10, 29, -19, 

-24, 43, -19, -31, -6, 37, -31, -34, 65, -31, -48) , 
and made the modified conjecture that this completes the list. 

The two other open cases are "dual" (see [AGKL]) to those two, and they would follow from them. 
Deconstructing the Ingenious Human Proof 

For the sake of definiteness and exposition, let's focus on the original Clark-Lewis [CL] recurrence, 

( x(n — l)-\-x(n) -r I i\ I / \ 

x{n + \) = \ if x(n-l) + x(n) is even; ^ 

\—x(n — V)-\-x{n), if x{n — V) + x(n) is odd. ' 

but following the proof strategy of [AGKL] . 

General Outline of the Amleh-Grove-Kent-Ladas Proof Strategy 

1. Try to prove that every trajectory is bounded. It would then follow immediately from the 
pigeon-hole principle, that there exists a pair of integers that shows up twice as {x{n),x{n + 1)), 
and once that happens it is trapped for ever into an orbit, but a priori infinitely many ultimate 
periods may show up, for different pairs of initial values. 
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2. In order to prove the boundedness assertion of 1 more specifically, look for two integers, ci and 
C2 such that, empirically for now: 

\x{n)\ < ci|a;(-l)| +C2|a;(0)| 
(for the present case ci = 1, C2 = 1 works). 

3. Try to prove {Fundlneq) by induction. You will probably fail at first, but whenever you fail, you 
can add more hypotheses, and try to prove a more general statement, proving several inequalities 
at once by induction, including the one you really need. You may have to keep adding more and 
more statements, and this process might never end (and indeed if (Funlneq) is false, it better not 
end!), but if in luck, this would end. We will describe it in more detail later. 

4. We now know that every trajectory must end-up with a cycle. We want to characterize all these 
cycles. Since a cycle is "cyclic", we can, without loss of generality, make the largest element (in 
absolute value) be the x{l) . Also since all the states are mod 2, we can take the initial conditions 
a;(— 1) = a and x{0) = b or x{—l) = a and x{0) = —b , where a and b denote (symbolic) positive 
integers. (We should separately treat the cases 6 = and a = 0, but this is really easy). 

By cyclicity (considering the last two elements x{n — l),x{n) as our a:;(— 1), x(0), and a;(l) as an 
entry that comes after), we have the necessary condition 

\x{l)\ < ci\x{n — 1)\ + C2\x{n)\ , 

that in addition to |x(n — 1)| < |x(l)|, \x{n)\ < |a;(l)|, enable us to "rule out" lots of possible parity 
sequences defined below. 

5. Define the parity sequence of a trajectory x{n) to be the sequence x{n) mod2. For example, 
the parity sequence of the trajectory 

11, 5, 8, 3, -5, -1, -3, -2, 1, 3, 2, -1, -3, -2, 1, 3, 2, -1, -3, . . . 

is 

1,1,0,1,1,1,1,0,1,1,0,1,1,0,1,1,0,1,1,... . 

We are looking for candidate parity sequences that do not violate the necessary condition of 4. 
Note that if we know the parity sequence of a trajectory, we can express each term of the trajectory, 
in particular the last two, as certain explicit linear combinations ofa = x(— l),6 = x(0). 

6. Using 4, we first empirically discover "all" the possible parity sequences of potential cycles. These 
parity sequences turn out (in all the successful cases encountered so far) to be regular expressions 
parametrized by few integer parameters. For example, for the present case they turn out to be 

(oii)™n™^ , (1011)1™^ , {oii)"'n"'^o , (oii)'"n'"^oi . 

7. Rigorously prove, by induction on mi,m2, ■ ■ ■, that these are the only ones, by using the 
inequalities of 4, but with symbolic mi,m2, ... (in addition to the symbolic a, b). 



(Fundlneq) 
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8. You can now express in closed form, in terms of a := x{—l) and b := x(0), the general two 
last terms, let's call them x{n — l),x{n), of such a trajectory. To investigate whether it can form 
a periodic orbit, try to solve the system of two linear equations 



x{n — 1) = a 



x(n) = b 



{orx{n — 1) = a 



x{n) = —b) 



getting a homogeneous set of two equations with two unknowns, a and b. Set the determinant of 
the 2x2 matrix of the resulting system of equations to zero, and get an expression in mi, m2, . . . 
that is usually (i.e. generically) not zero. For those rare cases where it is 0, you would get a hope 
for a cycle, and solving for a and b would give them to you. (This is a bit analogous to finding 
eigenvalues, followed by the corresponding eigenvectors). 

Before we can teach the computer how to discover and prove CoUatz-type theorems, we need to 
teach ourselves, in somewhat greater detail, the brilliant ideas of [AGKL] .Let's do this step-by-step 
for the original Clark-Lewis recurrence 



Proving Boundedness 

By using procedure FindPreScheme in the Maple package LADAS, the computer quickly makes the 
conjecture 



The natural approach would be to use induction on n. But it turns out that we (usually) need a 
stronger statement, so that we will have more elbow-room for the inductive argument. 

We will shortly describe how the computer can dynamically construct this stronger statement, but 
let us first describe how the computer (and of course, us) can rigorously prove its validity, once it 
has been formulated. We will give more details than in customary in human discourse, avoiding 
the phrases and words "it is easily seen", "obviously" (and whenever possible) "similarly". This 
level of detail is necessary if we want to be able to transfer our know-how to our silicon brethern. 

In the present case (recurrence (CL)), the stronger statement turns out to be as follows. Let 



A=\x{-1)\ + \x{0)\. 

Case I(n): If x{n — 1), x(n) are both odd then: 

\x{n - 1)1 < ^ , \x{n)\ < A . 
Case n(n): If x{n — 1) is odd and x{n) is even, then: 

\x{n-l)\<A , |x(n)|<^ , \ - x{n - 1) + x{n)\ < A 
Case Ill(n): If a;(n — 1) is even and ,x{n) is odd, then 

\x{n-l)\<A , |a:(n)|<^ , \ - x{n - 1) + x{n)\ < A 




if x{n — 1) -I- x{n) is even ; 
if x{n — 1) -I- x{n) is odd . 



{CL) 



x{n)\ < \x{-l)\ + \x{0)\ 



for n > —1 
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(In the present example of (CL), the two last cases can be combined into "if x{n — 1) + x{n) is 
odd" , since the three inequalities are identical. Usually this is not the case) . 

The natural approach would be to prove it by induction. There are three states 

/(n): (a;(n — l),x(ra))=(odd, odd) , 

II{n): {x{n — 1), a;(n))=(odd, even), 

III{n): (x(n — 1), a;(n))=(even,odd). 

Let's consider them each at their turn. 

If we are currently in state /(n) (i.e. (x(n — 1), x(n)) = (odd, odd)), then the next "state" is either 
I(n + 1) (i.e. (x{n),x{n + l))=(odd, odd)) or II{n + 1) (i.e. {x{n), x{n + l))=(odd, even)). 

Case la: (x(n — 1), 2;(n))=(odd, odd) , {x{n),x{n + l))=(odd, odd) . 

We have to prove statement I{n + 1), in other words: |x(n)| < A , \x{n + 1)| < ^4 . 

In terms of x{n — l),x(n) we have to prove: |a;(ra)| < A , ||a^(ra ~ 1) + ^3;(n)| < A . 

The first inequality is already contained in the (inductive) premise, while for the second we use the 
triangle inequality 

\^x{n - 1) + ^x{n)\ < ^\x{n - 1)| + ^|x(ra)| , 
By the inductive hypothesis we know that \x{n — 1)\ < A , \x{n)\ < A, so it follws that 

\\x{n-l) + ]^x{n)\ < \\x{n-l)\ + ^\x{n)\ <]^A + ]^A<A , 
since | + | < 1- 

Case lb: (a;(n — 1), x(n))=(odd, odd) , (x(n), x(n + l))=(odd, even). 
We have to prove //(n + 1), in other words: 

\x{n)\<A , \x{n+l)\<A , | - x(n) + x(n + 1)| < vl . 
Expressed in terms of x(n — 1) and x{n) these are 

\x{n)\<A , \{x{n-l)+x{n))/2\<A , \ - x{n) + {x{n - 1) + x{n))/2\ < A . 
Cleaning up, we have to prove: 

|a;(n)|<^ , \^x{n-l) + ^x{n)\<A , \^x{n - I) - ^x{n)\ < A . 



The first inequality is contained in the inductive hypothesis, the second one is identical to the one 
done above, while, for the third one, we have, once again by the triangle inequality: 

\^x{n ~ 1) - ^x{n)\ < \^\A + \-^\A < A , 

since |^| + |^| < 1. 

If we are currently in state II(n), the next state is necessarily III(n+l). We have to prove 
III(n+l): \x{n)\ < A, \x{n + 1)| < ^, | - x{n) + a:;(n + 1)| < A. 

In terms of x(n— 1), a;(n) these are: |x(n)| < A,\—x{n—l)+x(n)\ < A,\—x{n)—x{n—l)+x{n)\ < A. 

Simplifying, these are: \x{n)\ < vl, | — x{n — 1) + x{n)\ < A, \ — x{n — 1)| < A. 

The first two inequalities are part of the inductive hypothesis, while the third one follows from the 
deep fact that | — 1| = 1: 

\-x{n- 1)1 = |(-l)a;(n - 1)| = | - l||x(n - 1)| = l\x{n - 1)| = \x{n - 1)| < ^ . 

If we are currently in state Ill(n), the next state is necessarily II(n+l). In general, this would be 
different than the previous case, but in this example it is identical, so I hope that the reader will 
forgive me for using "similarly" , since this is not only similarly (in the colloquial human sense) but 
identical, and even a computer can realize that, and prevent a duplication of effort. 

How to teach the Computer to Argue as above 

In the above proof, we needed to know 

(i) If we are currently at a certain state{n), what state(s) is (are) next? 

(ii) Rewrite the expressions given in terms of x{n),x{n + 1), in terms of x{n — l),x(n), using the 
rules of the recursion, according to the current state. 

(iii) Know how to apply the triangle inequality, or realize that the inequality that we have to prove 
is already part of the inductive hypothesis. 

Tasks (i) and (ii) are obvious, while for (iii) we need to prove that a bunch of inequalities 

\Li\<A , \L2\<A , ... , \L^\<A 

imply another one 

|M|<A . 

If in luck, M is already one of the Lj's, or M = — Lj, for some i. Otherwise, we do a double do-loop 
looking for i and j such that one can write 

M = ciLi + C2Lj , 
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and then we check that |ci| + |c2| < 1. If this fails, we try a triple do-loop, looking ioi 1 < i < j < 
k < n such that 



|ci| + |c2| + |c3| < 1, etc. 

All this involves solving simple systems of linear equation, that Maple can do very fast. 

OK, A Computer can automatically prove, by Induction, such a "Scheme" of inequal- 
ities, but we sure need humans to come up with it! 

Wrong, of course! The computer starts with just what it really wants to know, namely \x{n)\ < A, 
and tries to prove it by induction. Chances are that there wouldn't be enough assumptions, so to 
paraphrase Guru Greg Chaitin, we "add it as an axiom", i.e. hypothesis, and keep trying to prove 
our "partial scheme" by induction, and whenever we need another assumption, let's just add it. If 
we are lucky, this process would eventually halt, and we would be done. 

Let's illustrate it with the (CL) recurrence. We start with the single inequality \x{n)\ < A in 
all states. Obviously wc need to add \x{n — 1)| < A right away. This suffices for the transition 
/(n) =^ I{n + 1), and /(n) =^ II{n + 1). Alas for II{n) =^ III{n + 1) and //I(n) =^ II{n + 1) we 
can't prove [ — x{n — 1) + x{n)\ < A, so we add it as part of the statement, and as we saw above, 
this suffices. 

In general, we may have to add more and more hypotheses, but if all goes well, we would eventually 
converge. 

Conjecturing the possible parity sequences of Cycles 

Procedure DAPf of our Maple package LADAS constructs all feasible parity sequences of a given 
rule, of any given length. This is a recursive procedure that starts out with all feasible symbolic 
trajectories of length 3, phrased in terms of a and b that conform to the convention that the third 
element, x{l), is the largest in absolute value, and using a and b to denote symbolic positive 
integers. For example for Rule 6 of [AGKL] 



the feasible trajectories (with the above convention that the third entry, x{l) is largest in absolute 
value) are: 



(here we use 2 instead of 0, since in the program denotes "either odd or even"). So we already 
know that it can't start with both entries being odd (since then x(l) wouldn't be able to be larger 
than both a and b). How does our computer generate the possible trajectories of length 4?. For 
each parity sequence of length 3 it considers the two extensions obtained by appending 2 (even) 
and l(odd). Then Maple automatically solves the system of linear inequalities (featuring absolute 



M = ciLi + C2Lj + csLfe 




x{n — 1) + x{n) is even ; 
x{n — 1) + x{n) is odd . 



{E6) 



[[a, -b, -a-b], [1, 2, 1]], [[a, -b, -a - b], [2, 1, 1]] 
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values, Maple can do it!) and if there are no solutions, discards that option. Step by step it 
constructs the potential parity sequences, and also expresses all the entries symbolically. 

After Maple gathered enough data, it looks for patterns, and by using procedure gGREplus (of the 
Maple package LADAS), the computer finds, empirically, all the regular expressions in the alphabet 
{0, 1} that are parity sequences of trajectories that do not violate any of the necessary conditions 
of 4. To prove that these are the only possibilities, rigorously, the computer first (symbolically) 
computes, (by using powers of matrices, or guessing, using the obvious format and then proving it by 
induction), explicit expressions for the last two elements of a trajectory x(n — l),x{n). For example, 
for (-B6) (resorting back to the convention of the paper, rather than the program, of denoting an 
even integer by rather than 2) there are three such possible symbolic regular expressions: 

(Oll)'^i , (011)"^U , (011)'"i01 . 

To prove that these are indeed all of them, simply use induction combined with the inequality of 4 
and rule out, for example the possibility (011)™iO, that would violate (for symbolic mi as well as 
symbolic a and b) that inequality. This too can be done automatically. 

Finding All the Cycles 

Now that we know what symbolic parity sequences may show up, we can explicitly express the last 
two terms of a trajectory (that is a potential cycle) in terms of a, b and mi,m2, ■ ■ ■■ Simply use 
symbolic powers of matrices (that Maple can easily do by finding the eigenvalues of the relevant 
matrices). For example, for rule (£'6) and the symbolic trajectory (011)™^, we have 

x{-l)=a , x{0) = -b , x{n - 1) = {-l)""' + (-1)"''+^ , x{n) = a^-^)""' . 

To investigate whether there exist mi (a positive integer) and positive integers a and b such that 
this forms a cycle, we have to solve: 

x{—l)=x{n — l) , x{0) = x{n) 

In other words 

a = a(-l)"^^ + (-l)'"i+^6 ,-b = a{~)"'' . 
Moving everything to the left: 

a(l-(-l)'"i) + (-l)"*i6 = , a{~)"''+b = . 

In order for there to be a non-zero solution (a, 6), the determinant: 

1 - (-l)'"^ + = , 

must vanish, but that can never happen. Similarly the other two (potential) symbolic trajectories 
can be ruled out. The only contributors are those where either a = or 6 = 0, and these easy cases 
are treated separately. 
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Output 



My beloved computer, Shalosh B. Ekhad, used my Maple package LADAS to consider, more generally, 
all 256 possible second-order difference equations of the form 

and x{n) is even ; 

and x(n) is odd ; 
and x{n) is even, 
and x(n) is odd. 

for (ai, . . . ,as) € {—1, 1}^. It was successful in 144 cases (more than a half!). All these computer- 
generated theorems (and proofs!), complete with the proving schemes for the fundamental inequal- 
ities (and their detailed proofs!), can be viewed from the webpage of this article 

http : //www ■ math . rutgers . edu/~zeilberg/mamarim/iiiamarimhtml/ collatz . html . 

The Maple package COLLATZ 

In addition to the main package LADAS, we also developed a purely empirical Maple package, 
COLLATZ, to conjecture the cycles of generalized Collatz-type transformations of the form n — > 
UiH + bi if n = i mod m, {i = 1 . . .m) such that mai and iai + bi are integers. 

Further Work 

This is but the tip of an iceberg, except that the specific weight of our ice is much closer to 1 than 
that of real ice. Obvious extensions would be to consider general parity-dependent transformations 
on Z^, rather than those of the form (a, b) (6, F{a, b)) considered here. More generally, we should 
consider transformations of the form 

ia,b)^F{a,b) , 

where F{a,b) = Fij(a,b) ii a = i mod m, b = j mod m, and we have linear transformations 
that map pairs (a, b) such that if a mod m and b mod m equal i and j respectively, then Fij(a, b) 
outputs a pair of integers. 

Why stop at transformations of Z^? We can also consider transformations of Z'^ for d = 3, 4, . . .. 

The two remaining open cases of [AGKL], that we can't do either, may very well be amenable to the 
present approach, except that the ci,C2 that feature in the inequality |x(n)| < ci|x(— 1)| -|-C2|x(0)| 
are rather large, and our program was not fast enough to conjecture them? It would be a good 
idea to take our amateurish and slow program and speed it up (say using Java and C), and perhaps 
be able to tackle these. On the other hand, we may need yet another idea, and perhaps the two 
remaining cases of [AGKL] are just as hard as the original Collatz problem. 

What about affine-linear transformations? The original Collatz recurrence is a very simple first- 
order recurrence given in that way. What about using the [AGKL] method, as computerized here. 



x(n + 1) 



/" aix{n — l)-\-a2x{nj 
2 ' 
a3a:(n — l)+CK4x(n) 
2 ' 

a5x{n — 1) -|- aQx{n) 
, aTx{n — 1) -|- asx{n) 



if x{n — 1) 

if x(n — 1) 

if x{n — 1) 

if x(n — 1) 



IS even 

is odd 

is odd 

is even 
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to systematically explore even one-dimensional analogs of the 3x + 1 problem, given by higher 
moduli, rather than two? See the Maple package COLLATZ mentioned above, for very preliminary 
(empirical) investigations. 

And who knows? Perhaps a sufficiently higher-order recurrence, amenable to the present approach 
(or to a yet-to-be-discovered extension) would imply, via an appropriate specialization, the good-old 
Collatz problem? 

So there is a lot to do. Research is a relay race. Myself (and even my beloved Shalosh) are already 
tired of this project, but we do hope that other people (and machines!) will take over. Amen. 
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